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Abstract 

In this work, some new integrable and nonintegrable cosmological models of the Hofava- 
Lifshitz gravity are proposed. For some of them, exact solutions are presented. Then these 
results extend for the F(R) Hofava-Lifshitz gravity theory case. In particular, several inte- 
grable cosmological models of this modified gravity theory were constructed in the explicit 
form. 

1 Introduction 

More than one year ago Hof ava proposed a theory, the so-called Hofava-Lifshitz quantum gravity, 
which is a power-counting renormalizable theory with consistent ultra-violet behavior pQ. In this 
theory, the scaling of the gravitational system at short distances exhibits a strong anisotropy 
between space and time 

x l ^bx\ t^b z t. (1.1) 

This important relation between space and time coordinates can be realized with the Arnowitt- 
Deser-Misner decomposition of the metric of the form 

ds 2 = -L 2 dt 2 + g ij (dx i + L l dt)(dx j + L 3 dt), (1.2) 

where gij is the spatial metric (roman letters indicate spatial indices), L and Li are the lapse 
and shift functions, respectively. Recently the F(R) Hofava-Lifshitz quantum gravity has been 
proposed [2] (see also [3]-[8]). In this work, we consider integrable aspects of the usual and F(R) 
modified Hofava-Lifshitz gravity theories ( see e.g. [T0]-[TT]h 

The paper is organized as follows. In section 2, we study some integrable and nonintegrable 
FRW cosmological models of Hofava-Lifshitz gravity. In the next section 3, we extend some of 
these results for the modified F(R) Hofava-Lifshitz quantum gravity case and present its some 
exact integrable models. Exact solutions of some integrable models of the usual Hofava-Lifshitz 
gravity were considered in the section 4. Section 5 is devoted to the conclusion. 

2 Cosmological models of Hofava-Lifshitz gravity 

In this work, we restrict ourselves to the detailed balance case. In this case, the action of Hofava- 
Lifshitz gravity takes the form (z = 3): 

/o 2 2 ijk 

dttfxJgNl-siKyKV - \K 2 ) + ^ ( - ^—R a V J R L k + 
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%ff J - 0/0 7 iA (^-^ 2 +Afl-3A 2 )], (2.1) 



where Vj are the covariant derivatives defined with respect to the spatial metric gij, t % ^ k is the 
totally antisymmetric unit tensor and K,X,p,w = consts. Here K^j and C lJ are the extrinsic 
curvature and the Cotton tensor, respectively, which are given by 

Kij = O.bL- 1 ^ ~ ViLj - VjLi), C lJ = g-^^VkiBi ~ 0.25BS}). (2.2) 

As our main interest in this work is the cosmological aspects of HL gravity, we impose the pro- 
jectability condition. Consider FRW spacetime with the scale factor a(t) and metric 



ds 2 = -dt 2 +a 2 (t) 

that is 



dr 2 



1 — kr 2 



(2.3) 



L = l, g l} : = a 2 (t) 7ii , L l = 0, (2.4) 

where k can take any value but it is related to (— , 0, +) curvatures according to sign. In this case, 
the Friedmann equations we can write 
in the H-form 

p = -a(2H + 3H 2 ) +/i(o), (2.5) 
p = 3aH 2 + f 2 (a), (2.6) 
p = -3H(p + p) (2.7) 

or in the N-form 

p = -a{2N + 3N 2 ) + f 1 {a), (2.8) 
p = 3aN 2 + f 2 {a), (2.9) 
p = -3N{p + p). (2.10) 

Here N = In a, H = N and 

fi = A) + h^ 2 + faa-\ f2=m+V2a- 2 +ma-\ (2.11) 
where a = 2k~ 2 (3A — 1) and 

00 = 0.75^ 2 A 2 a- 1 , f3 2 = -0.5kp 2 Aa~\ f3± = -0.25p 2 k 2 a-\ (2.12) 

?7 = -0.75p 2 A 2 a- 1 = — j9o, m = l.bkp 2 ka- 1 = -3/? 2 , m = 0.75/x 2 fc 2 'a^ 1 = -3/3a. (2.13) 

In the next sections we will study some integrable aspects of the Friedmann equations in the H-form 
(2.5)-(2.7) or in the N-form (2.8)-(2.10) and try solve some of them. 

2.1 Integrable models 

Let us start from the presentation of some examples of integrable HL cosmological models. We 
assume that N, a, H satisfy one of Painleve equations so that we obtain 30 new integrable HL 
models. Consider examples [below a, (3,^,5, k and p are arbitrary constants] (see also [T0]-|11|). 

1) P/ - models. In our cosmological case, this model as and other P-type models have 5 
particular submodels that means 5 type cosmological models. For this case we have: 
i) Pia - model. The Pja - model we write in the following closed form 

p = -a{2H + 'iH 2 ) + h{a), (2.14) 

p = 3aH 2 + f 2 {a) 7 (2.15) 

N = 6N 2 +t, (2.16) 

p = -3H(p + p) (2.17) 
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or 



p = -a(2N + 3N 2 )+f 1 (a) 7 (2.18) 

p = 3aN 2 +f 2 (a), (2.19) 

N = 6N 2 +t, (2.20) 

p = -3N(p+p). (2.21) 

We believe that the system (2.14)-(2.17) [or its equivalent (2.18)-(2.21)] is integrable. Note that 
the case a = 1, fi — corresponds to General Relativity (GR). Keeping in mind that the Pia - 
model has the form (2.14)-(2.17) [or equivalently (2.18)-(2.21)], for short, we here write it as 



N = 6N 2 +t. (2.22) 



We also note that in the systems (2.14)-(2.17) or (2.18)-(2.21), the equation (2.22) plays the role 
of the EoS p — p(p). Finally we would like to note that similarly the other models we below write 
in the same short form as (2.22). 

11 1 U - - . W-l /"^ /-I S~^l I * 



ii) P/B - model: 
hi) Pic - model: 

iv) Pjd - model: 

v) Pie - model: 

2) P a - models. 

i) Piia - model: 

ii) Pub - model: 

iii) Puc - model: 

iv) Pud - model: 

v) Piie - model: 

3) Pj/j - models. 

i) Piiia - model: 

ii) Piub - model: 

iii) Pine - model: 

iv) Phid - model: 

v) Piiie ~ model: 



h = (Sa 2 +t. (2.23) 

H = 6H 2 + t. (2.24) 

H aa = 6H 2 + a. (2.25) 

H NN = m 2 + N. (2.26) 

N = 2N 3 +tN + v. (2.27) 

a = 2a 3 + ta + v. (2.28) 

H = 2H 3 + tH + v. (2.29) 

H aa = 2H 3 + aH + v. (2.30) 

H NN = 2H 3 + NH + v. (2.31) 

11 A 

N = —N 2 --(N-aN 2 -/3)+~fN 3 + —. (2.32) 

11 s 

a= -a 2 - -(a- aa 2 - (3) + 7a 3 + -. (2.33) 
at a 

H = ±H 2 -j(H-aH 2 -f3)+~fH 3 + ^. (2.34) 

H aa = \-H 2 a --{H a - aH 2 -P)+ ~fH 3 + A. 
Ha tl 

H NN = ^H 2 N - ±(H n - aH 2 -fi) + 1 H 3 + A 



(2.35) 
(2.36) 
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4) Pjy - models. 

i) Piva - model: 



ii) Pivb - model: 

iii) P/yc - model: 

iv) Pivd - model: 

v) Pive - model: 



1 ^ 

N = — N 2 + 1.5N 3 + UN 2 + 2(t 2 - a)N + — . (2.37) 

1 s 

o = —a 2 + 1.5a 3 + 4ta 2 + 2(t 2 - a)a + -. (2.38) 

1 ^ 

ff = — H 2 + 1.5H 3 + AtH 2 + 2(t 2 - a)H + — . (2.39) 

H aa = t^H 2 + 1.5H 3 + 4aH 2 + 2(a 2 - a)H + A ( 2 .40) 



1 s 

H NN = —H 2 N + 1.5H 3 + 4NH 2 + 2(N 2 - a)H + —. (2.41) 
2H ti 



5) Py - models. 

i) Pva - model: 



= (^ + J^jW 2 -- t (N- 7 N) + t- 2 (N l) 2 (aN + + 6N jf_\ ^ ■ (2.42) 

ii) Pvb - model: 

a = (^ + -^-)d 2 - I( a - 7 a) + i" 2 (a - l) 2 (aa + /3a" 1 ) + 6 < a+1 \ (2 .43) 

iii) Pyc - model: 

H =(^ + J^j)H 2 -\{H- 1 H)+ r 2 (H l) 2 (aH + (1H^) + ^L±12 . (2 . 4 4) 

iv) Py£> - model: 

Haa = (^ + ^)ff 2 - ±(2T a - 7#) + a" 2 ^ - if '(off + /JIT" 1 ) + SH ^_\ 1) . (2.45) 

v ) Pve - model: 

Hnn = (^ + Jj^)H% - jj(H N - 7 H) + N- 2 (H \) 2 (aH + /Sif" 1 ) + S _W_Lll, (2 . 46 ) 

6) Py/ - models. 

i) Py/A - model: 



N N-l N-tJ \t t- 1 N-t, 

+ t~ 2 (t - 1)- 2 N(N - 1)(N - t) [a + ptN- 2 + 7 (t - 1)(JV - 1)~ 2 + tfi(t - 1)(JV - t)~ 2 ] . (2.47) 
ii) Py/s - model: 

a = 0.5 ( - + -^-r + — — I a 2 - f - + — !— ■ + 



K a a— 1 a — t) \t t—1 a — t, 
+ t~ 2 (t - ir 2 a(a - l)(a - i) [a + ^a" 2 + 7 (t - l)(a - 1)~ 2 + - l)(a - t)~ 2 ] . (2.48) 



iii) Py/c - model: 

/ I) 5 ( 

i? H-tJ \t t-1 H-t 



^^o.5[i +T A +T AU 2 -ri + -L +T7 i-i^ 
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+ r 2 {t - iy 2 H{H -1)(H- t) [a + ptH- 2 + 7(f - 1)(H - l)~ 2 + St(t - - t)~ 2 ] . (2.49) 
iv) Pvid - model: 

H aa =0.5(1 + + J—\ H 2 -( 1 - + ^ + J—\ H a 



H H-l H-aJ a \a a-l H - a / 
2 {a-l)- 2 H{H-l){H-a) [a + f3aH- 2 + 7(0 - 1)(H - l) -2 + Sa{a - 1)(H - a)' 2 ] . (2.50) 



+ a 

v) Pvie - model 



1 1 1 \ „o / 1 1 1 



; : I77 ' 77— r ' 77— vj /y v l.- ' t-7 ' 77-v ) ' 

+ iV" 2 (iV - iy 2 H(H -1)(H- N) [a + [3NH- 2 + j(N -1)(H- l)" 2 + 5N(N - 1)(H - N)- 2 ] . 

(2-51) 

Note that all Painleve equations can be represented as Hamiltonian systems that is as (see e.g. 
jlOj and references therein) 

i = f , (2.52) 
f . -f q , (2,3, 

where E(q,r,t) is the (non-autonomous) Hamiltonian function. Consider some examples (see e.g. 
|10] and references therein). 

1) P/-models. In this case, q,r,F read as 

q = r, (2.54) 
r = 6q 2 +t, (2.55) 
E = 0.5r 2 - 2q 3 - tq. (2.56) 

2) P/7-models. In this case we have 

q = r-q 2 -0M, (2.57) 
r = 2qr + a + 0.5, (2.58) 
E = 0.5r 2 - (q 2 +0M)r- (a + 0.5)q. (2.59) 

3) P///-models. In this case we have 

tq = 2q 2 r - k 2 tq 2 - (29 1 + l)q + kit, (2.60) 
tf = -2qr 2 + 2k 2 tqr + (29 1 + l)r - fc 2 ((6» 1 + 8 2 )t, (2.61) 
tE = q 2 r 2 - [k 2 tq 2 + {2e l + l)q-k l t]r + k 2 {e l +e 2 )tq. (2.62) 

So in this subsection we presented 24 new HL cosmological models. Note that as integrable 
systems, these models admit n-soliton solutions, infinite number commuting integrals of motion, 
Lax representations etc. 



2.2 Nonintegrable models 

Here we consider some known and new HL models induced by some ODEs. These ODEs are 
nonintegrable so that the corresponding HL cosmological models are nonintegrable (see also |llj). 

1) ACDM cosmology. We start with the ACDM cosmology. Here we present 5 submodels 
[We remark that in fact only the Ai - model corresponds to ACDM cosmology] . 

i) Ai - model: 

N = 0.5A- 1.5iV 2 . (2.63) 

ii) A 2 - model: 

a = 0.5A- 1.5a 2 . (2.64) 
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iii) A 3 - model: 

H = 0.5A- 1.5H 2 . (2.65) 

iv) A 4 - model: 

Haa = 0.5A-1.5Hl. (2.66) 

v) A 5 - model: 

H NN = 0.5A- l.hH%. (2.67) 

2) Pinney cosmology. It induced by the Pinney equation. Let us present 5 submodels, 
i) Pinneyi - model: 



where £ = £(t). k = const. 
ii) Pinney2 - model: 



iii) Pinney3 - model: 

iv) Pinney4 - model: 

v) Pinney 5 - model: 



iv=£(i)7V + A (2.68) 



k 

'd = £(t)a+—. (2.69) 



H = Z(t)H+± (2.70) 



H aa =t(a)H+^. (2.71) 



H NN =C(N)H + ^. (2.72) 

3) Schrodinger cosmology. For this model also write 5 submodels. 

i) Schrodinger i - model: 

N = uN + kN, (2.73) 

where u = u[t), k = const. 

ii) Schrodinger - model: 



iii) Schrodinger 3 - model: 

iv) Schrddingcr4 - model: 

v) Schrodingers - model: 



d = ua + ka. (2.74) 
H = uH + kH. (2.75) 
H aa = uH + kH. (2.76) 

H NN = uH + kH. (2.77) 

4) Hypergeometric cosmology. Let us we present 5 submodels. 

i) Hi - model: 

N = t-^l - t) _1 {[(« + + l)t- "f]N + apN}. (2.78) 

ii) H 2 - model: 

a = - i) _1 {[(o! + /? + l)i - i\a + afia}. (2.79) 

iii) H 3 - model: 

H = t~\\ - i) _1 {[(Q! + + l)t- ~f]H + apH}. (2.80) 
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iv) H4 - model: 

H aa = a-\l - a)- x {[(a + /3 + l)a- j]H a + apH}. (2.81) 

v) H5 - model: 

Hnn = N- X {1 - N)- 1 ^ + /3 + l)N- j}H N + a(3H}. (2.82) 

Finally we would like to note that the above presented HL models, in particular, have some 
solutions which describe the accelerated expansion of the universe. 

3 Cosmological models of F(R) Hofava-Lifshitz gravity 

Let us now consider the F(R) Hofava-Lifshitz gravity. Following [2], its action we write as 

S = 0.5k- 2 J dtd 3 Xy /gNF(R), (3.1) 

where k 2 = 16irG and 

R = K l} K ij - \K 2 + R + 2^iV ' M (n M V ' v n v - n u V^) - L (3) 5y . (3.2) 

The case F(R) = R corresponds to the original Hofava-Lifshitz gravity. In this section, we con- 
struct integrable and nonintegrable F(R) Hofava-Lifshitz gravity models induced by some ODSs. 
For the metric 

ds 2 = -L 2 dt 2 + a 2 (t){dx 2 + dy 2 + dz 2 ), (3.3) 

the scalar R takes the form 

R = 3(1 - 3A + 6^)H 2 L- 2 + 6^£T MlnO^" 1 )]*- (3-4) 
In the FRW case and as L = 1, the equations of motion read as [2] 

= F — 2(1 — 3A + 3/i)(iJ + 3H 2 )f' — 2(1 — 3\)RF" + 1[iR 2 F" + RF" ) + K 2 p m , (3.5) 
= F -6[{l-3\ + 3[i)H 2 + [iH]F' +6[iHRF" -K 2 p m , (3.6) 

where p m and p m are the pressure and energy density of a perfect fluid that fills the Universe. 
3.1 Integrable models 

Some classes integrable models can be constructed with the help of known integrable ODEs. It is 
the main idea of the work (see also 11 ). One of famous representatives of such integrable ODEs 
are Painleve equations. It is the case that we are going to use to get integrable Hofava-Lifshitz 
gravity models. Let us demonstrate it. The lazy way to do it is the assumption that the function 
F(R) or its sister f(R) — F(R) — R satisfy some integrable ODEs, in our case, one of Painleve 
equations. Note that there are 6 integrable Painleve equations with each of them can be relate 

4 integrable F(R) Hofava-Lifshitz models so that totally we have 30 models. Now we are in the 
position to present these equations (or models) [F = dF/dR, F = dF/dt etc]. 

1) Fi(R) - models. 

a) Fja{R) - model: 

/" 

or 

F" 

b) F IB (R) - model: 

} 

or 



6f 2 + R 


(3.7) 


6F 2 + R. 


(3.8) 


6f 2 + t 


(3.9) 


6F 2 + 1. 


(3.10) 
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c) F IC (R) - model: 



or 



d) F ID (R) - model: 



or 



e) F IE (R) - model: 



or 



2) Fn(R) - model, 
a) Fha{R) - model: 



or 



b) Fu B (R) - models: 



or 



c) Fjjc{R) - models: 



or 



d) F IID (R) - models: 



or 



e) Fjje{R) - models: 



or 



3) Fin {R) - models. 

a) Fiha(R) - model: 



or 



faa = 6f + a 

F aa = 6F 2 + a. 
Inn = 6/ 2 + N 



F NN = 6F 2 + N. 



Shh = 6f + H 



F HH = 6F 2 + H. 



f =2f + Rf + a 
F" = 2F 3 + RF + a. 
f = 2f +tf + a 
F = 2F 3 + tF + a. 
faa = 2f + af + a 
F NN = 2F 3 + NF + a. 

f NN = 2f + Nf + a 



F NN = 2F 3 + NF + a. 



Jhh = 2f + Hf + a 



F HH = 2F 3 + HF + a. 



f" = - f f' 2 -j(f'-af-P)+ 1 f+° 1 
F" = ±F' 2 -l(F'-aF 2 ~P) +1 F 3 + l 
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b) Fhib(R) - model: 



or 



c) F UI c{R) - model: 



or 



d) F IIID (R) - model: 



or 



e) F IIIE {R) - model: 



or 



4) F IV (R) - models, 
a) Fjva(R) - model: 



or 



b) F IVB {R) - model: 



or 



c) F IVC {R) - model: 



or 



d) F IVD {R) - model: 



or 



c) Five(R) - model: 



F=^-\{F-aF*-P)+ 1 F* + ^. 



faa = J. f a \{fa af - /3) + 1 f + j 

F aa = ^F 2 a - 1 -(F a -aF 2 -P)+^ + ^. 



Inn = jf% - ±(f N - af - P) + 7 / 3 + j 
F NN = jF 2 - 1(F - aF 2 - /3) + 7 F 3 + |. 



Ihh = j.f 2 H ^{Ih af 2 - /3) + 7 / 3 + j 
F HH = ±F 2 -± (F a - aF 2 f3) + 7 F 3 + | . 



/" = ^/' 2 + 1-5/ 3 + 4Rf 2 + 2(R 2 - a)f + - f 

1 A 

F" = —F 2 + 1.5F 3 + 4RF 2 + 2(R 2 - a)F + — 

2_r r 



f=^f 2 + 1.5f + 4tf + 2(t 2 -a)f + ^ 

1 A 

F = —F 2 + 1.5F 3 + 4tF 2 + 2(t 2 - a)F + -. 

Zr t 



faa = Yfti + L5 / 3 + 4 «/ 2 + 2 (« 2 - «)/ + J 

Faa = ^F 2 a + 1.5F 3 + AaF 2 + 2(a 2 - a)F + |. 



Fnn = ^F 2 + 1.5.F 3 + ANF 2 + 2{N 2 - a)F + - 



.fen = ^/Ir + 1-5/ 3 + AHf 2 + 2(H 2 a)f + - f 



or 

Fhh = ^pF% + 1.5F 3 + AHF 2 + 2(H 2 - a)F + |. (3.46) 

5) F V (R) - models. 

a) F V a(R) - model: 

/" = (-j + j^j)/ 2 - 7/) + R- 2 (,f l) 9 (a/ + + (3.47) 

or 

F " = { h + t-i )f ' 2 i {F ' 7F) + k ~ 2{F l)2{aF + + ^T- + i 1} - (3 - 48) 

b) F VB (R) - model: 

/ = + J^)/ 2 - \{f if) + t- 2 (f l) a (a/ + M- 1 ) + SJ ^ 1 (3-49) 

or 

F = { h + F^i )p2 ¥ F - lF) + r2{F 1)2{aF + pF ~ 1] + 6F F F - + i 1} - (3 - 50) 

c) F VC (R) - model: 



faa = + J^y)./ 2 - kfa if) + a" 2 (/ - 1) >/ + /J/" 1 ) + 



(3.51) 



or 



F aa = (-L + ^R 2 - " 7*0 + a" 2 (F - 1) 2 (^ + Z^" 1 ) + (3-52) 
d) F VD (R) - model: 

/jviv = (^ + j^j)f 2 N ±Vn If) + a- 2 (.f l) 2 (af + Pf- 1 ) + Sfi f f ^ i 1) (3.53) 

or 

F NN = (-L + - ^ - 7^) + N- 2 (F l) 2 (aF + PF^) + ^Llll . (3.54) 

c) Fve{R) - model: 

fHH = (-j + y-^)/! - ^(.fe - 7/) + H- 2 (.f l) 2 (af + Z?/" 1 ) + (3 ' 55) 

or 

Fhh = + - 7^ - 7^) + ^" 2 (^ - 1) 2 (^ + Z^" 1 ) + (3.56) 

6) - models. 

a ) Fvia(R) - model: 

+*-»(*-i)-V(/-i)(/-» lo + wr' + iW - W - I)"' + *«(* - l;!/ •«> " . «■■"••> 

Instead of this equation we can consider the following one 



f /-I f-RJ' \R R-l f-R 



F F-l F-RJ \R R-l F-R 
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+ Br 2 (R - 1)- 2 F{F - 1)(F -R) a + PRF' 2 + j{R - 1)(F - 1)~ 2 + SR(R - 1){F - R)~ 

(3.58) 

b) F VIB {R) - model: 



J /-i f-tJ V t-i f-t / 

+ t- 2 (t - l)- 2 f(f - 1)(/ - t) [a + (3tf- 2 + 7 (i - 1)(/ - l)- 2 + St(t - 1)(/ - ty 2 ] . (3.59) 

+ r 2 {t - l)- 2 F{F - l)(F -t)[a + /3tF- 2 + j(t - l)(F - 1)~ 2 + 5t(t - l)(F - ty 2 ] . (3.60) 
c) Fvic{R) ~ model: 

1 1 1 \ .Q /l 1 1 



faa = 0.5 - + — — + fa-[- + 7 + ~ f fa 

\f /-I f-aj \a a-1 f-aj 

+ a- 2 (a - l)- 2 f(f - 1)(/ - a) [a + paf' 2 + 7 (a - 1)(/ - 1)~ 2 + Sa(a - l)(f - a)" 2 ] . (3.61) 
and 



F F-l F-a/ a V a a- 1 F ~ a / 
+ a" 2 (a - 1)" 2 F(F - 1)(F - a) [a + /3aF~ 2 + 7(0 - l)(F - 1)~ 2 + 5a(a - l)(F - a)' 2 ] . (3.62) 
d) F VID (R) - model: 

fNN = °' 5 (7 + J^T + 7^v) ^ " (at + aT^i + jhi) fN 

+ N- 2 (N - l)- 2 /(/ - 1)(/ - N) [a + /37V/- 2 + 7 (^V - l)(.f - 1)~ 2 + 6N(N - 1)(/ - TV)" 2 ] . 

(3.63) 

and 

1 1 1 \ o / 1 1 1 



+ N~ 2 (N - 1)- 2 F(F - 1)(F -N)[at + /WVF~ 2 + j(N - 1)(F - 1)~ 2 + SN(N - 1)(F - TV)" 2 ] . 

(3.64) 

e) F V ie{R) - model: 

+ H- 2 (H - iy 2 f(f - 1)(/ - if) [a + l3H.r 2 + 1 (H - 1)(/ - 1)~ 2 + 6H(H - l)(.f - H)' 2 ] . 

(3.65) 

and 

' 1 



f »« = - 5 {f + f—i + f^) f »-{h + —i + f-h) f » 

+ H- 2 (H - l)- 2 F(F - l)(F -H)[a + f3HF~ 2 + 7 (if - 1)(F - 1)~ 2 + SH(H - l)(F - H)- 2 ] . 

(3.66) 

So we presented new 30 HL cosmological Fj(R) models (J — I, II, III, IV, V, VI) which are 
integrable due to integrability of Painleve equations. 

Exact solutions of integrable HL models. All above constructed integrable HL models 
admit (may be infinity number) exact solutions. Let us here present some of them [for simplicity, 
we give just some particular solutions for the sister function f(R) = F(R) — R and only for some 
models] (see e.g. [TU]V 

1) The Fha(R) - model. 
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1) As our first example, we consider the Fjja(R) - model (3.17). Let us present its some 
solutions. For example, this model has the following particular solutions: 

f(R) = /(^;1.5) = V-(2V 2 +^)- 1 , (3.67) 

f(R) ee f(R;l) = -l, (3.68) 
H 

1 3 /? 2 

f(R) = f{R;2) = -- w - l , (3.69) 

f(R) ee f(R-3)= 3k " QR2 (& + 1Q ) f3 70) 
i? 3 + 4 fl B + 20i? 3 -80' 

f(i?) ee f(&4) = 1 I 6 ^ 2 ^ 3 + 1 °) 9i? 5 (i? 3 + 40) 

V ' ' R i? 6 + 20# 3 -80 # 9 + 60^ + 11200' 

f(R) = /(^;0.5e) = -eV (3.72) 

and so on. Here 

i) = (In 4>) & , <P(R) = C\Ai(-2- 1/3 R) + C 2 Bi(-2- 1/3 R), (3.73) 

Ci = consts and Ai(x),Bi(x) are Airy functions. 

ii) Similarly, for the Fub(R) - model (3.19) we have the following particular solutions 

f(R) = /(t;1.5)=^-(2^ 2 +t)- 1 , (3.74) 

f(R) = 1, (3.75) 

/(i?) ee /(t ;2 ) = I-^ i) (3.76) 

f{R) = /(t ' 3) -^T4"t6 + 20i3-80' (3J7) 

= f(fi)- 1 I 6t2 ft 3 + 1 °) 9f 5 (t 3 +40) 

/W ~ ' ; ^ 4 ^ t + te + 20i 3 -80 + 60t 6 + 11200' ( ' 

f(R) = /(i;0.5e) = -e^ (3.79) 

and so on. Here 

iP = (ln^) t , 0(t) = dAi^-^h) + C 2 Bi(-2- 1 'H). (3.80) 

2) The F IIIA (R) - model. 

i) Our next example is the Fjha(R) - model (3.27). It has the following particular solutions: 

f(R) = f{R;v 1 ,0,0,-vn4) = v 2 y/k, (3.81) 

f(R) = f(R; 0,-2^,0, 4^ - 4) = R[\n{R^f + In (R»*e^)}, (3.82) 

f(R) = f(R;-v*v 2 ,0,v 2 1 (vl-v 3 v i ),0)= — , (3.83) 

V3R 2 " 1 + v 2 R Vl + va 

f(R) = f(R; 2v x + 3, -2v x + 1, 1, -1) = - + ^ - , (3.84) 

J? + V\ + 1 

/(-R) = f{R;v u -v 1 i4,v 3 ,-v 3 V2)=V2, (3.85) 

/(i?) = / = - £l (ln^ (3.86) 

and so on. Here 

ip(R) = R v [CiJ v { y /e^R) + C 2 Y u {^eT^R% {e it C l = consts), (3.87) 
and J v {x), Y v {x) are Bessel functions. 
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ii) The corresponding particular solutions of the Fhib{R) - model (3.29) are given by 

f(R) = f(t;v 1 ,0,0,-v 1 vl) = v 2 yi, (3.88) 



f{R) = /(t;0, -21^,0, 4 I /iiA 2 -i/|)=t[ln(^) 2 + In (t^e^)], (3.89) 



f(R) ee f(t; -v\v 2 , 0, v\{vl - ^ 4 ), 0) = , (3.90) 

/(i?) = /ft; 2^ + 3, -2v x + 1, 1,-1) = t + Vl , (3.91) 

< + i/i + 1 

/ft) = f{t\vi,-vivl,vz,-vavt)=v*, (3-92) 
/(£) = / = - ei (ln^) t (3.93) 

and so on. Here 

<p(t) = + C 2 Y u {y/wt)]. (3.94) 

Similarly we can present the exact solutions of the other models. 

3.2 Nonintegrable models 

Let us here present some known and new HL models induced by some ODEs. These ODEs are 
nonintegrable so that the corresponding HL cosmological models are nonintegrable. Consider 
examples. 

i) We start from models induced by the hypergeometric differential equation. We write here 5 
versions of this model. 

R-version: 

R(l - R)f" + [c - (a + b + l)R]f - abf = 0. (3.95) 

It has the solution f(R) — 2 F(a, b; c; R) which is the hypergeometric function, 
t-version: 

t(l - t)f + [c - {a + b + l)t]f -abf = (3.96) 

with the solution f(t) = 2 F(a,b;c;t). 
a-version: 

o(l - a)f aa + [c - (a + b + l)t]f a - abf = (3.97) 

with the solution f(a) = 2 F(a 1 b;c;a). 
N-version: 

N(l - N)f NN + [c - (a + b + l)t]f N -abf = (3.98) 

with the solution f(N) = 2 F(a, b; c; N). 
H-version: 

if (1 - H)f HH + [c - {a + b + l)t]f H -abf = (3.99) 

with the solution f(H) = 2 F(a, b; c; H). 

ii) Another example is the case when f(R) satisfies the Pinney equation 



/ +eiW + ^V=0. (3.100) 



t2_(R) 

P 

If £i = 1, ^2 = K — const, this equation has the following solution (see e.g. [11] ) 

f(R) — cos 2 R + k 2 sin 2 R. (3.101) 
Its "t-form" is f(t) — cos 2 t + k 2 sin 2 t which is the solution of the Pinney equation 

/ + &(*)/+ ^T=0 (3.102) 

as £i = 1, £2 = k = const. 

iii) Let us we present one more example. Let the function / satisfies the equation 

/" = 6/ 2 - 0.5 52 , (g 2 = const) (3.103) 
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or 

/ = 6/ 2 - 0.5g 2 . (3.104) 
These equations admit the following solutions 

f(R) = p(R) (3.105) 

and 

f(R) = f(t) = p(t), (3.106) 

where p(R) and p(t) are the Weierstrass elliptic functions. 

Finally we would like to note that the above presented HL models, in particular, have some 
solutions which describe the accelerated expansion of the universe. 



4 Cosmological solutions 

It is important find exact solutions of HL gravity models (see e.g. [I]- [5]). In our case, all above 
presented HL models admit exact solutions. It is important that some of these solutions describe 
accelerated expansion of the universe. Let us present some cosmological solutions of some above 
presented HL models. As an example, consider the P/j-models that is the equations (2.27)-(2.31). 
These equations have the following particular solutions (see e.g. |10]-|11]). 
i) P IIA - model (2.27): 

N(t) = N(t;v 1 ,Q,Q,-v 1 i>$)=i>2\fi, (4.1 
N(t) = N(t; 0,-2^,0, 4^i ia> - v\) = <[(ln<^) 2 + m (t^e" 2 )], (4.2 

N(t) = N(t; -v\v 2 , 0, - 0) - ^ - - , (4-3 

N{t) = N(t; 2v! + 3, -1v x + 1, 1, -1) = -j + V \ (4.4 

t + V\ + 1 

N{t) = N(t;u 1 ,-y 1 vlu 3 ,-u 3 4)=v 2 , (4.5 

N(t) = N = -aQxLtp)^ (4.6 

Hence we get the corresponding expressions for the scale factor a(t). We have 

a(t) = a(t; u u 0,0, -v x v\) = e v * r \ (4.7 
a(t) = o(t; 0,-2i^,0,4^ -v\) =t V3t e U2t e<- latVi;rt f, (4.8 



a{t) = a{t- 1 -v(v 2l Q,vi{v^-v 3 ,v i ) 1 Q)^e"^ 2 " 1 +^" 1 +^ 1 (4.9 

a(t) = <z(f; 2^i + 3, -2z/i + 1, l,-l) = e I +^ T , (4.10 

a(t) = a(t; 1/1,-^1/2^3,-1/31/2) = e" 2 , (4.11 

a(t) = a = e- £l(lnv)t . (4.12 



ii) Pub - model (2.28): 



a(t) ee a(t;i/i,0,0,-Uiui)=i/ 2 ^i, (4.13 

a(t) ee a(t; 0,-2^,0, 4i/ x i*j - z/f) = <[ln(<^") 2 + In (t" 3 e V2 )], (4.14 

^1—1 

a(f) = a(t; -i/fe, 0, v\{v\ - 1*1/4), 0) = , (4.15 

o(t) = aft; 2^ + 3, -2^ + 1, 1,-1) = f + ^ , (4.16 

t + Vl+l 

a(t) ee a{t;v x ,-viv\,V3,-V3v\) = v 2 , (4.17 

a(t) ee a = -ei(hup)t. (4.18 
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iii) P//c - model (2.29): 



H(t) 
H(t) 



H(t; 0,-2^,0,41/11/2 - v\) = t[ln {t^f + In (t^e" 2 )}, 



(4.19) 
(4.20) 



H(t) = H 
H{t) = H 



H(t) 



H(t) 




(4.21) 



(4.22) 



(4.23) 
(4.24) 



and so on. These exact solutions correspond to the different cosmologies. Some of them describe 
the accelerated and decelarated phases of the Universe. Finally we note that similarly we can 
present exact solutions of the other HL models constructed in the previous sections. 

5 Conclusion 

In this work, a new class integrable and nonintegrable cosmological models of the Hofava-Lifshitz 
gravity were proposed. For some of them, exact solutions are presented. To construct integrable 
models we use the well-known integrable systems, namely, Painlevc equations. Then these re- 
sults extend for the F(R) Hofava-Lifshitz gravity theory case. In particular, several integrable 
cosmological models of this modified gravity theory were constructed in the explicit form. 
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